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The central arm spectrometers for the PHENIX experiment at the 
Relativistic Heavy Ion Colhder have been designed for the opti- 
mization of particle identification in relativistic heavy ion collisions. 
The spectrometers present a challenging environment for event re- 
construction due to a very high track multiplicity in a complicated, 
focusing, magnetic field. In order to meet this challenge, nine dis- 
tinct detector types are integrated for charged particle tracking, 
momentum reconstruction, and particle identification. The tech- 
niques which have been developed for the task of event reconstruc- 
tion are described. 
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1 Introduction 



The PHENIX experiment [1] at the Relativistic Heavy Ion Collider (RHIC) 
is designed to measure hadrons, leptons, and photons produced in nucleus- 
nucleus, proton-nucleus, and proton-proton collisions at beam energies of up 
to 100 GeV/A with the primary goal of detecting the Quark-Gluon Plasma 
(QGP) and characterizing its physical properties. PHENIX consists of four 
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Fig. 1. Layout of the PHENIX detector showing ah four spectrometer arms within 
the PHENIX experimental hall. The collision point is in the center of the figure. 



spectrometers, or arms (see Fig. 1), including two muon arms designed pri- 
marily for muon identification at forward rapidities, and two central arms 
focusing on hadron, electron, and photon identification near mid-rapidity (see 
Fig. 2). This document describes the software used for pattern recognition, 
momentum reconstruction, and particle identification for the two central arm 
spectrometers of the PHENIX detector, which began taking data in May 2000 
until September 2000, a period which will be referred to here as Run 2000. 

This document is organized as follows. Section 2 provides a general overview 
of the detectors in the central arm spectrometers. Sections 3 through 10 dis- 
cuss the software for each individual detector component. Section 11 explains 
the momentum reconstruction and track model definition method. Section 
12 explains the inter-detector association algorithm. Section 13 explains the 
technique for particle identification using the time-of-flight detector. Section 
14 summarizes what has been discussed. 
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PHENIX Detector - First Year Physics Run 




□ Installed 
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Fig. 2. Schematic diagram of the PHENIX Central Arm Spectrometers. The abbre- 
viations are defined in the text. The darker shading indicates which detectors were 
implemented for Run 2000. 

2 The PHENIX Central Arm Spectrometers 

In order to provide a spectrometer with optimized particle identification ca- 
pabilities for various particle species, the central arm spectrometers, which 
each cover ±0.35 in pseudorapidity (j]) and span 90° in azimuthal angle (0), 
necessarily contain many components. When there are difi^erences in the two 
spectrometers, they will be referred to as the east and west arms. The posi- 
tions of the detectors will be described in a coordinate system with the origin 
at the center of the central arm magnet with the z-axis along the beam. This 
section will briefly describe each component moving outwards radially from 
the beam pipe. 

The innermost detector in each spectrometer is the drift chamber (DC) [2]. 
The drift chamber provides the primary vector and momentum measurement 
for charged particles traversing the spectrometers. The drift chamber consists 
of 40 planes of wires arranged in 80 drift cells placed cylindrically symmetric 
about the beam fine. Each drift chamber spans 90° in 0, has a radial sensitive 
region from 2.02 to 2.46 meters, and covers -80 < z < -|-80 cm. The wire planes 
are placed in an X-U-V configuration in the following order (moving outward 
radially): 12 X planes (XI), 4 U planes (Ul), 4 V planes (VI), 12 X planes (X2), 
4 U planes (U2), and 4 V planes (V2). The U and V planes are tilted by a small 
±5° stereo angle to allow for full three-dimensional track reconstruction. The 
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field wire design is such that the electron drift to each sense wire is only from 
one side, thus removing most left-right ambiguities everywhere except within 2 
mm of the sense wire. The wires are divided electrically in the middle at z = 0. 
The occupancy for a central RHIC Au+Au collision is about two hits per 
wire. The magnetic field bends particles in the x-y plane. Generally, the drift 
chamber active volume lies outside of the axial magnetic field, however there 
is a residual field within the volume that provides less than 1*^ of deficction 
in the bend plane. There is no bending, to first order, in the non-bend plane, 
except for tracks very close to the magnet pole edges or at low momentum [3]. 
The deflection of the tracks is sufficiently low that the track model assumption 
can be straight line trajectories in both planes. 

Attached to the back of the drift chamber, at an inner inscribed radius of 248 
cm, is a single layer of pad chamber called PCI. This is a pixel-based detector 
[4,5] with effective readout sizes of 8.45 mm in z by 8.40 mm in the x — y plane, 
providing three-dimensional space-point information for charged particles in 
the spectrometer. 

Behind PCI is a cylindrically shaped Ring Imaging Cherenkov Detector (RICH) 
[6,7]. The RICH is the primary detector for electron identification in PHENIX. 
It is a threshold Cherenkov detector with a high angular segmentation to cope 
with the high particle density expected in the most violent collisions at RHIC. 
There are two identical RICH detectors, one in each PHENIX arm. Each RICH 
detector consists of a gas vessel, thin refiection mirrors, and photon detectors 
made up of arrays of photo- multipher tubes (PMT's). During Run 2000, CO2 
was used as the Cherenkov radiator. Pions below 4.9 CcV/c do not produce a 
signal in the RICH, while electrons above 18 MeV/c emit Cherenkov light in 
the gas radiator. The Cherenkov photons are refiected by the mirrors in the 
RICH, and are focused onto the photon detectors. A charged particle track 
is identified as an electron when a sufficient number of RICH PMT hits are 
associated with it. 

On the west arm only, a second pixel pad chamber, named PC2, is located 
behind the RICH at an inner inscribed radius of 419 cm. This detector was 
not installed for Run 2000. PC2 has an effective readout size of 14.25 mm in 
z by 13.553 mm in the x — y plane. 

On the east arm only, a Time Expansion Chamber (TEC) [8] is placed behind 
the RICH with a radial active area from 4.1 to 5.0 meters. One of the four 
installed TEC detector planes were instrumented for Run 2000 providing az- 
imuthal coverage oi (f) — 45°. The TEC contains four planes of wire readout 
(expandable to six planes in the future). The wires are divided electrically at 
their center at z=0. The readout electronics consist of a 5-bit non-linear fiash 
ADC. The TEC provides 2-dimensional tracking in the x-y plane, contributes 
to electron- pion discrimination via dE/dx measurements, and provides high 
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Fig. 3. The TOF detector system mounted on the PHENIX East arm in the PHENIX 
experimental hall showing 10 panels of the detector. 

momentum resolution for high pt particles. The TEC is located outside of the 
axial magnetic field, so a straight hne track model assumption is valid for TEC 
track reconstruction. 

Placed at an inner inscribed radius of 490 cm on each arm is a third layer of 
pixel pad chamber, called PC3, with an effective readout size of 16.7 mm in 
z by 16.0 mm in the x — y plane. Only the east arm of PC3 was installed for 
Run 2000. 

Located on the east arm only is the time-of-flight detector (TOF) placed 
at an inner inscribed radius of 510 cm (between PC3 and the calorimeter) 
and spanning 30° in 0. The TOF detector, shown mounted on the east arm 
in Fig. 3, serves as the primary particle identification device for charged 
hadrons in PHENIX. The TOF wall is finely segmented in order to cope with 
the expected highest multiplicity events in central Au + Au collisions, which 
is about 1000 charged particles per unit rapidity. One TOF panel contains 
96 plastic scintillation counters (Bicron BC404) with photomultipher tubes 
(Hamamatsu R3478S) at both ends of the slat. A total of 10 TOF panels, 960 
slats, and 1920 channels of PMT's were installed for Run 2000. Each slat is 
oriented along the r-0 direction and provides timing as well as longitudinal 
position information of particles that hit the slat. 

On the east arm only, placed behind the TOF detector, are two planar arrays 
of Lead-Glass calorimeters (PbGl) imported from the CERN SPS experiment 
WA98 [9], as were the TOF detectors. Behind PC3 on the remainder of the 
east arm he two planar arrays of Lead- Scintillator calorimeters (PbSc) [10]. 
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Both of these arrays span 45° in (j) on the east arm. Covering the full 90° 
in (p on the west arm. and placed behind PCS, are four PbSc arrays. The 
calorimeters, collectively referred to as EMC, provide photon identification 
information, particle energy, and time-of-flight measurements. 

Trigger, timing, and collision centrality information is provided by a set of 
detectors including a Beam-Beam Counter (BBC) [11], a pair of Zero- Degree 
Calorimeters (ZDC) [12,13], and a Multiplicity Vertex Detector (MVD) [14]. 
The BBC consists of two identical sets of counters installed on both sides of 
the interaction point along the beam. Each counter is composed of 64 PMT's 
equipped with a quartz Cherenkov radiator and covers the 3.0 < rj < 3.9 region 
with full azimuthal coverage. The ZDC's are small transverse area hadron 
calorimeters which measure the energy of unbound neutrons in small forward 
cones {9 < 2 mr) around each beam axis. One ZDC is located on either side 
of the interaction region and consists of 3 modules, each with a depth of 2 
hadronic interaction lengths and read out by a single PMT. Both time and 
amplitude are digitized for each of the 3 PMT's as well as an analog sum of 
the PMT's for each ZDC. 

The MVD is a highly segmented sihcon strip and pad detector. It is designed to 
measure and trigger on the number of charged particles, as well as to measure 
the collision vertex with high accuracy. The detector has a large pseudorapid- 
ity coverage (approximately -2.64 < rj < 2.64) and full azimuthal coverage, 
making it possible to study charged particle production as a function of 77 
and (f) on an event-by-event basis. The detector is composed of two parts: two 
concentric barrels, or shells, spanning -32 < 2; < 32 cm, and two endcaps at 
z=±35 cm. Each barrel shell consists of six rows of detectors, where each row 
is formed by 12 panels of silicon strip detectors with a 200 iim pitch. The 
disk-shaped endcap sections are each composed of six silicon pad wafers each 
segmented into 252 pads. When fully instrumented, about 35,000 channels will 
be read out from the MVD. 



3 Beam-Beam Counters and Zero Degree Calorimeters 

The Beam-Beam Counters (BBC) and Zero Degree Calorimeters (ZDC) are 
designed to define the start timing of the nuclear interaction for time-of-fiight 
measurements, to provide the trigger signal of the nuclear interaction, and to 
provide a coarse measurement of the collision vertex along the beam axis. The 
following equations are used to calculate the start timing (Tq) and the vertex 
position along the beam axis {Z^ertex) for both detectors: Tq = (Ti + T2)/2 
and ^uertea; = (^1 — ^2)/2c, where Ti and T2 are the average hit timing for 
each set of BBC counters, and c is the velocity of light. 
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Fig. 4. The difference in the measured BBC and ZDC Zyertex positions as a function 
of the measured BBC Zygrtex position from the Run 2000 data. 

The average hit timing in each BBC set is defined by taking a truncated 
average of the timing over the PMT's that have a hit (defined as a PMT with 
a vahd timing signal) out of the 64 possible PMT's. In the case that there 
is no peak structure in the hit timing distribution in either BBC set due to 
low particle multiplicity, the arrival times of the first leading particle from the 
beam collisions are used as the hit timing value. 

The ZDC pulse height is converted to hadronic energy equivalent using pe- 
ripheral interactions with low neutron multiplicity in the ZDC where single- 
and double-neutron peaks are clearly seen. The fitted single-neutron peak po- 
sition is equivalent to 65 AGeV per beam for Run 2000. This fit procedure 
was performed using an analog sum channel since pedestal offset corrections 
are reduced in this channel. The offset arises because the analog-to-digital 
converters used for the BBC and ZDC are self-gating, so a minimum signal 
must be present to start the digitization. 

The time-of-ffight to the ZDC was calculated using the PMT digitized times 
after correcting for pulse height slewing via a look-up table. The typical signal 
rise-time is 2 nsec. The time digitizer least-count scale was calibrated with 
standard delay cables. The resulting time measurements arc used to calculate 
To and Zi,crtex as outlined above using a C++ class whose data members are 
stored for each event for both the BBC and ZDC. Fig. 4 plots the difference 
between the measured ZDC and BBC vertex positions as a function of the BBC 
vertex position. The resolution in the difference is 2.0 cm, which is dominated 
by the contribution from the ZDC time resolution. 
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4 Multiplicity Vertex Detector 



The Multiplicity Vertex Detector (MVD) is designed to measure and trigger on 
the number of charged particles produced in high multiplicity events in Au+Au 
collisions and to measure the collision vertex with high accuracy. Vertex finding 
is performed by the MVD barrel section where the hit information from the two 
shells of detectors can be combined and projected onto the z-axis. Multiplicity 
measurements are performed by both the barrel and pad sections. For a vertex 
at the center of the MVD (z — 0), the inner shell covers the range -2.55 < rj < 
2.55 and the Si pad detectors cover the range 1.79 < I77I < 2.64 at each end 
of the detector. Both the inner shell and the pads have complete azimuthal 
coverage. 

The MVD event processing, written in C-|— 1-, is handled by an MVD event 
reconstruction object, which is created and evaluated for each event by the 
main MVD reconstruction module. The event reconstruction object contains 
data members holding the hit information for each event and the final vertex 
and dN/dr) results. The content of these data members are evaluated by a 
call to the internal event processing method. Strips, pads, hit clusters, and 
the vertex are treated as individual objects for the given event. The content 
of the cluster objects are generated from the strip objects. Strip, pad, and 
cluster objects are accessed by a vertex object class which holds the methods 
containing the main reconstruction algorithms. 

Two different algorithms are used to find the vertex position. The first algo- 
rithm is denoted the pseudo-tracking method. This vertex algorithm takes each 
pair of clusters (where a cluster is a group of adjacent channels in the same 
row, shell, and panel of the MVD which all have signals above a threshold) 
in the inner and outer shells and treats every combination in the same row of 
the MVD barrel as a potential track. This combination of track candidates is 
projected back to the beam axis to determine the vertex expected for a track 
which hit these two locations. For each of these possible vertex locations, a 
channel in a histogram is incremented. Most pairs of hits do not correspond 
to a real track. However, when all pairs are considered, the true vertex loca- 
tion appears more frequently than the other locations. A clear peak appears 
in the histogram corresponding to the vertex position, as shown in Fig. 5 
from simulated central Au+Au events using the HIJING event generator [15]. 
When this algorithm finds the correct vertex location, it is correct within a 
standard deviation, a, of about 190 /xm for simulated central events and 224 
//m for simulated minimum bias events. At high occupancy (around 45% in 
the inner shell) this method becomes more inefficient (less than 98%) and a 
different approach must be used. 

The alternative vertex algorithm, used for high multiplicity events, is the cor- 
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Fig. 5. MVD vertex finding histogram for the pseudo-tracking method using simu- 
lated central HIJING Au-|-Au events. The narrow peak defines the vertex. 

relation method. Unlike the pseudo-tracking method, the correlation method 
makes use of the ADC information from the individual channels. The algo- 
rithm is based upon the assumption that the pattern of energy loss per chan- 
nel in the inner and outer shells of the MVD arc similar. The pattern of ADC 
values in the outer shell of the MVD is scaled down by the ratio of the inner- 
to- outer shell radius and stored in an array. The scaled pattern yields an offset 
proportional to the vertex location when compared to the unsealed inner shell 
pattern. By finding the value of the offset which makes the two patterns most 
similar, the vertex position is found. A histogram which plots the difference 
between the pattern of hits in the inner and outer shells will show a sharp 
minimum which corresponds to the location of the vertex, as shown in Fig. 
6 from central HIJING Au-|-Au events. This method gives the correct vertex 
with a standard deviation of 233 /im for simulated central events and 250 /im 
for simulated minimum bias events with an efficiency greater than 99% for 
inner shell occupancies above 5%. 

The algorithms for calculating dN/ drj are different in the MVD barrel section 
and endcap sections. In the MVD barrel, only the inner shell is used. In cen- 
tral Au+Au collisions, the occupancy in the inner shell is high [14], therefore 
the algorithm makes no attempt to find individual hits. Instead, the barrel is 
divided into groups of 64 adjacent strips which are in the same row (azimuthal 
segment). For each of these groups of channels, a dN/dr] value is calculated 
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Fig. 6. MVD vertex finding histogram for the correlation method using simulated 
central HIJING Au+Au events. The narrow dip defines the vertex. 

as follows: the total ADC value for all channels in this group is added up, the 
range of pseudorapidity (A77) is calculated, the pseudorapidity at the center 
of the group of strips is calculated, a geometric correction (for non-normal 
incidence) is applied, the number of particles is estimated by dividing the 
corrected ADC signal by the ADC signal expected for a minimum ionizing 
particle (mip), and finally dN/df] is estimated as the number of particles di- 
vided by Arj. The average dN/dr] for a set of events is calculated by taking 
the average in each bin of rj. 

In order to calculate dN/dr] from the MVD pads, the algorithm assumes that 
a single particle hits only one pad. The occupancy of the pad detectors is 
around 15-20% for central events. The algorithm estimates the number of 
particles associated with the observed ADC value in an individual pad using 
an ADC distribution which is taken from the distribution corresponding to 
the ADC response expected for a single particle at normal incident angle in 
low multiplicity events, and the occupancy of the pad detector. The algorithm 
does not associate an integer number of particles with each hit. Instead a mean 
number of particles which would be associated with a given ADC channel is 
calculated. This number can then be converted into a dN/ dr] value using the 
vertex location and the pad geometry. 
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Fig. 7. The pad chamber pixel pattern. The three pixels form a cell at the center of 
the figure. Nine connected pixels form a pad. 

5 Pad Chambers 

A new pad geometry for cathode readout of MWPC's was designed for the 
PHENIX pad chambers [5,16]. The basic component in this geometry is a 
pad consisting of nine smaller connected copper electrodes, called pixels (Fig. 
7). The pad is read out by a single preamplifier and discriminator. The pads 
are interleaved in a repeated pattern. The area defined by three adjacent 
pixels belonging to three different pads is called a cell. This arrangement 
saves a factor of three in the number of readout channels with marginal loss 
of performance. 

Pad chamber information, which is stored in the PHENIX Objectivity database, 
can be separated into three different groups: geometry, high voltage and elec- 
tronics. All groups are accessed via object-oriented classes in the offline anal- 
ysis, using time stamps as the primary key. The geometry group holds infor- 
mation on the geometry parameters of the chambers themselves, e.g. wire- 
spacing, as well as survey information, while the high voltage group keeps 
track of when trips have occurred. The electronics group contains threshold 
and gain calibrations, as well as information on where hot or inactive channels 
are located. 

The list of fired pads from an event is read by the hit reconstruction algorithm, 
which first translates it into a list of fired cells. The limit on the number of 
allowed inactive or hot pads in a cell is an input parameter whose default 
value is one, i.e. at least two out of three pixels in a cell have to be part of 
alive channels. From the cell list, cluster objects are built for each chamber by 
identifying all neighboring cells. The size of each cluster is checked against a 
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Fig. 8. The distribution of the number of cells per cluster in PCS from a) simulation 
and b) data. 

pre-defined table to estimate the number of traversing particles that created 
the cluster. The default settings of this table were determined from simula- 
tions [16]. If the number of particles in a cluster is considered to be greater 
than one, the cluster will be recursively split depending on its shape until only 
one-particle clusters remain. Once a cluster created by a single particle has 
been extracted, the average hit position is calculated by averaging over the 
coordinates of the cells in the cluster. The entire procedure is repeated until 
the hst of fired cells is empty. 

In order to determine the reconstruction efficiency, a pad chamber response 
simulator, which simulates which pads fired in an event, is utilized. Given the 
coordinates of a particle entering and exiting the pad chamber, the location of 
the avalanche along the closest wire is calculated. The particle may also give 
rise to additional avalanches on neighbouring wires. The total avalanche charge 
is obtained from a Landau distribution. The induced charge distribution on 
the pad cathode plane is then calculated from an empirical formula [17,18]. 
Finally, when all particles of the event are processed, noise charge drawn from 
a normal distribution is added to each pad in the chamber. If the total collected 
charge on a pad exceeds the threshold value, the pad is registered as fired. The 
thresholds used in the simulation were set to match the cluster sizes (Fig. 8) 
and the very high efficiency of the real chambers. 

The position resolutions, presented in Table 1, are determined by comparing 
the reconstructed hit positions with the positions where the original GEANT [19] 
particles passed through the chambers. The position resolutions for PCI are 
shown in Fig. 9. By simulating minimum ionizing particle pairs with small 
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Fig. 9. The simulated position resolution along the wire (left panel) and across the 
wire (right panel) in PCI. 
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Fig. 10. The probability to reconstruct two simulated particles as two separate hits 
as a function of the separation of the particles in PCI both along the wire (left 
panel) and across the wire (right panel). 

separation between them, the probability for resolving the nearby tracks was 
calculated. Shown in Fig. 10 is the two-track resolution along and perpendicu- 
lar to the wires for the PCI chamber. The resolution is defined as the distance 
where 50% of the tracks are separated by the reconstruction. The system can 
be optimized for two-track resolution at the expense of an increased number 
of ghost hits (due to falsely split clusters) or loss of single track efficiency by 
operating at lower sensitivity (leading to reduced cluster sizes). The study 
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Table 1 

Pad chamber position and two-track resolutions. 



Simulation Results 


PCI 


PC2 


PCS 


Position resolution along wire [mm] 


1.6 


2.7 


3.2 


Position resolution across wire [mm] 


2.3 


4.0 


4.8 


Two-track resolution along wire [cm] 


2.4 


4.0 


5.0 


Two-track resolution across wire [cm] 


2.9 


4.6 


5.3 



presented here results in a negligible number of ghosts. While the detection 
efficiency as measured with cosmic rays is greater than 99.5%, the reconstruc- 
tion efficiency is slightly lower as determined by projecting reconstructed drift 
chamber tracks from minimum bias data into PCI. This results in a recon- 
struction efficiency calculated to be better than 98%, as expected from the 
pad chamber simulation. 



6 Drift Chamber 



Track reconstruction within the drift chamber is performed using a combina- 
torial Hough transform (CHT) technique [20,21]. In this technique, the drift 
chamber hits are mapped pair-wise into a feature space defined by the polar 
angle at the intersection of the track with a reference radius near the mid- 
point of the drift chamber, 0, and the inclination of the track at that point, 
a. The a variable is proportional to the inverse of the transverse momentum, 
thus facilitating limited searches for specific momentum ranges and providing 
an initial guess for the momentum reconstruction procedure. Fig. 11 provides 
a schematic illustration of these variables. 

This coordinate mapping results in a gain effect that greatly enhances the 
peak-to-background ratio in the feature space. If n is the number of points 
on the track, the nominal peak height is n{n — l)/2. In order to minimize 
computation time, a and were only calculated for hit pairs separated in 
azimuthal angle by physically reasonable amounts corresponding to pt — 150 
MeV/c, comparable to the natural resolution of the spectrometer. Fig. 12 
shows an example of a region of the drift chamber hits and the associated 
feature space from a simulated HIJING central RHIC Au-|-Au collision. In 
order to reconstruct tracks in regions where some wires have been removed for 
internal supports, the algorithm first looks for tracks that traverse both the 
XI and X2 wire regions, then it looks for the remaining tracks that traverse 
the XI or X2 regions only. 

After reconstructing a track in the magnetic field bend plane, the direction of 
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Fig. 11. Illustration of the Hough transform parameters for drift chamber track 
reconstruction. The outline shows the drift chamber active volume. The circles rep- 
resent drift chamber hits along the particle trajectory. 




Fig. 12. The left panel shows simulated hits from a central Au+Au collision for 
a small physical region of the drift chamber. The right panel shows the Hough 
transform feature space for this region. Tracks appear as peaks in this plot. 

the track is specified by and a. Track reconstruction in the non-bend plane 
is first attempted by integrating information from PCI reconstructed clusters, 
which contain z information with a resolution of 1.89 mm. A straight line 
projection of the bend-plane track is made to the PCI detector and a road 
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Fig. 13. The reconstruction efficiency in HIJING-simulated central Au+Au events 
as a function of the transverse momentum of the track. 

is defined about that projection. If there is an unambiguous PCI association, 
the non-bend vector is defined by the PCI cluster and the z-coordinate of 
the event vertex. This assumption is not vahd for secondary tracks, but those 
cases are later examined during momentum reconstruction. If there is no PCI 
cluster association, or if there are multiple PCI association solutions, the non- 
bend vector is determined from the stereo wires. For this procedure, a plane, 
called the X-plane, is defined by (p and a, but extending in z. The path of the 
track trajectory through the X-plane is a straight line and is determined by 
including the information measured by the stereo wires. To do this, all drift 
chamber hits are defined as lines in physical space parallel to the anode wire 
and perpendicular to the drift direction. All stereo wire hit lines are intersected 
with the X-plane with the intersections belonging to the track defining a line 
in the X-plane. Again, a Hough transform is used to extract these intersections 
by using the z coordinate at which the track crosses the reference radius, zed, 
and the polar angle at that point, 0. Since the X-plane is defined by a single 
track, a unique solution in this space is extracted. 

The method of determining the drift chamber track reconstruction efficiency 
will be discussed in more detail in Section 12. Using this method, the drift 
chamber tracking efficiency as a function of transverse momentum is shown in 
Fig. 13 with less than 1% spurious (ghost) tracks. The resulting spatial reso- 
lution of the tracks is comparable to the single-hit resolution of the detector. 
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7 Time Expansion Chamber 



Track reconstruction in the Time Expansion Chamber (TEC) is performed 
using a Combinatorial Hough Transform method using the identical variables, 
a and 0, and with an identical method for track reconstruction in the magnetic 
field bend plane as the drift chamber track reconstruction. 

A detailed simulation of the TEC detector response is used to determine the 
TEC tracking efficiency. This simulation calculates the average number of 
primary ion collisions for each track section that traverses a single drift cell 
[22] . Fluctuations of the number of primary collisions are modelled by Poisson 
statistics [23]. Secondary ionization is calculated by generating a spectrum 
of 5 electrons [23]. The total number of primary and secondary ion pairs are 
calculated and the distance of the electron cluster to the wire, along with the 
drift time, are calculated assuming a straight line drift. Gas gain fluctuations 
are simulated using a 'polya distribution [24] . After the addition of noise, and 
the addition of the electronics shaping time response, the charge collected on 
the wire as a function of time bin is calculated. The charge is finally converted 
into FADC channels for each track. The final FADC contents are obtained by 
summing the contributions from all charged tracks in the event. The track re- 
construction efficiency is determined by embedding individual simulated tracks 
processed through the response simulation into a real data event, reconstruct- 
ing the resulting merged event and determining how many of the embedded 
tracks are correctly reconstructed. The reconstruction efficiency is shown as a 
function of the charged track multiplicity in the TEC in the left panel of Fig. 
14. The spatial resolution of reconstructed tracks in the TEC, shown in the 
right panel of Fig. 14, is determined directly from data taken in Run 2000 by 
comparing the residuals of tracks reconstructed independently in alternating 
planes of the TEC and fitting them to a Gaussian curve. The resulting resolu- 
tion, also confirmed within the simulation procedure, is determined to be 380 



8 Electromagnetic Calorimeters 

Pattern recognition in both the PbSc and PbGl Electromagnetic Calorimeters 
(EMC) includes single electromagnetic shower identification including calcu- 
lation of the shower position, energy, time-of- flight, overlapping shower sepa- 
ration, and hadronic shower rejection. The technique is based upon matching 
energy- and impact-angle-dependent descriptions of the projected electromag- 
netic shower shape as measured in a test beam of electrons in the momentum 
range 0.3-5.0 GeV/c [26] and reproduced in a detailed GEANT simulation of 
the detector response. An illustration of the shape of a shower is shown in 
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Fig. 14. The TEC track reconstruction efficiency as a function of the TEC charged 
track multipHcity obtained by evaluating simulated tracks embedded in real data 
events (left panel), and the spatial resolution of 380 /xm of the TEC based upon 
data taken in Run 2000 (right panel). The arrows on the left panel correspond to 
the 5% most central, and the 60-80% peripheral events in the Au+Au collisions. 



Fig. 15 for a particle impact orthogonal to, and at a 20*^ impact angle to, the 
calorimeter surface. 

The shower reconstruction algorithm includes two steps. First, a cluster search 
is performed. The cluster is defined as a contiguous group of calorimeter towers 
with measured energies above a certain threshold. At this point, the clusters 
are considered to be totally independent. Second, the clusters are split into 
subclusters associated with local maxima. A local maximum is defined as a 
tower which has a deposited energy in excess of the energies measured in any 
of the 8 surrounding towers. A subcluster (peak region) is formed from the 
towers of the cluster sharing the energy in each tower among photons located 
in each peak region. Here, an iterative procedure based on a parameterized 
shower profile is used [25] . Each subcluster is assigned a value for shower 
identification purposes. 

The simplest algorithm for determining the particle position in a cellular 
calorimeter is to measure the center-of-gravity of the shower, Xcg'- 

Xca ^ (1) 



where Xi is the x-coordinate of the center of tower i, and Ei is the energy 
deposited in the tower. For electron and photon position reconstruction, a 
parametrization of the shower center-of-gravity in units of tower width is ap- 
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Fig. 15. Energy deposited in a PbSc calorimeter tower as a function of the distance 
between the tower center and the impact point (left panels), the tower center and the 
shower center-of-gravity (right panels); a),b),d) and e) are for orthogonal impact; 
c) and f) are for a 20^^ impact angle. 
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plied [28]: 

^^^^2 sinh{l/2h) '"'i^Mmp." "1^2 



XcG = ^ \ ' ^ + ^; \X^m, + A - (51 < ^. (2) 



where is the particle impact point, h is the shower cross-sectional width, 
A is the mean displacement of the calculated shower center-of-gravity from 
the impact point Xi^p, and (5 is a phase shift related to the skewed shape of 
the shower projection. All of these parameters ( 6, A, and 5) are parameterized 
as a function of the incident particle energy, and angle of incidence, 9. A 
is given by 

A = Le// • sin{9), 



where Lg// is the effective depth of the shower penetration in the calorimeter 
as determined by the position of the cascade-curve median in the longitudinal 
direction (oc ln{E)) [29]. The value of b is given by 

b^bo + bi(E)sin'^(e), 



where bo is a constant and bi{E) oc ln{E). 5 is found to be essentially energy- 
independent and can be parameterized as a function of the impact angle only. 
Finally, the position of the particle impact on the calorimeter surface, Ximp, 
is calculated according to the formula obtained by inverting Equation (2): 

Ximp — b ■ arcsinh[2 ■ {Xcg — ■ sinh{l/2b)] — A + S, 
\XcG-S\<^. 

The calorimeter position resolution is well described by [27,28] : 

ax{E, 9) = ax{E, 0°) © {d ■ sin{9)), 



where (Tx{E, 0°) is the position resolution for an orthogonal incidence. The sec- 
ond term in this parametrization describes the contribution of the longitudinal 
shower fluctuations on the projection to the calorimeter surface. The value of 
d is O.SXo, where Xq is the radiation length in the calorimeter module. Using 
the electron test beam, a position resolution of h.lmm/ ^ E{GeV) -\- 1.6mm 
was achieved. 

A number of factors are taken into account when reconstructing the photon 
and electron energy in the calorimeters: the longitudinal energy leakage, the 
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light attenuation in the calorimeter module, the response dependence on the 
particle impact position, and the shower overlaps in the high multiplicity en- 
vironment, some of which will be discussed here. The nonlinearity of photon 
and electron energy measurements in the calorimeters is primarily connected 
with light attenuation in a calorimeter module. This is well described by the 

expression Eemc/E ^ E^att^ where Eemc is the energy measured in the 
calorimeter, and E is the actual photon or electron energy in GeV. 

In a high multiplicity environment, shower overlaps can seriously affect the 
electromagnetic shower energy measurements, thus degrading the energy res- 
olution and inducing systematic shifts in the measurement. For a tower energy 
threshold of 3 MeV, the mean cluster size in the PbSc calorimeter initiated 
by a 0.5 GeV (4 GeV) photon is 7 (19) towers. The larger the cluster size, 
the higher the probability for an overlap. However, the shower size can be 
reduced to 4-5 towers in the shower core. Towers belonging to the shower core 
are defined using a parameterized shower profile: 



where r is the contribution of the tower to the total shower energy i?, and 
El'^'^'^ is the predicted energy deposited in the tower. The e = 0.02 param- 
eter is chosen in order to have as few towers in a cluster as possible while 
keeping the energy resolution on the same level. The shower core defined in 
this manner contains about 92% of the total shower energy, dropping to 90% 
for a 20*^ impact angle. Using the shower core for energy measurements intro- 
duces an additional constant term of about 2% in the energy resolution, which 
slightly worsens the energy measurements of isolated photons and electrons, 
but helps to considerably improve the energy measurements in the high mul- 
tiplicity environment, where an energy resolution for the PbSc calorimeter of 



^ E{GeV) © 1.9% was achieved in electron test beams at the BNL AGS 
and the GERN SPS [26]. 

The time-of-fiight (TOP) information provided by the calorimeters (using the 
BBG as the start counter) is a useful tool for particle identification. This helps 
by enhancing the e/vr separation at low momenta (< 0.8 GcV/c), where the 
cncrgy-to-momentum comparison, as well as the shower profile analysis, lose 
effectiveness for c/tt discrimination. In addition, this is the primary tool for 
identifying neutral hadrons, which are a major source of background in the 
photon spectrum in the energy region of approximately 1 to 3 GeV. Also, the 
magnet poles and structural elements of the detectors in front of the calorime- 
ter are sources of non-vertex photons, part of which (~ 15%) can be rejected 
using the TOP information. Pinally, in the case of high multiplicity events, re- 
quiring consistency of the TOP values measured in the different towers within 
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Fig. 16. The distribution of hit PMT's around the projected ray from a reconstructed 
track in zero-field runs after the geometry cahbration. 

the same cluster can help to identify accidental shower overlaps. Several cor- 
rections are applied to measure the TOF, including slewing corrections and 
flight path corrections based upon the collision vertex, after which a resolution 
of better than 700 psec for the PbSc calorimeter was achieved. 



9 Ring Imaging Cherenkov Detector 

Electron identification using the Ring Imaging Cherenkov Detector (RICH) 
is performed with two C++ classes. The first class handles the geometry of 
the RICH detector by modelling the position, shape, and orientation of the 96 
mirror panels and the position of the 5120 PMT's in the RICH detector. This 
class also handles all optical tracing calculations. The second class uses the 
geometry class and performs the electron identification. The electron identi- 
fication class accepts a previously reconstructed track and searches for RICH 
PMT's that can be associated with the track. This is done by first refiecting 
the track about the RICH mirror as if it were a light ray, and then projecting 
the reflected track onto the RICH PMT array. If the track is an electron, the 
reficctcd ray should be projected to the center of a ring formed by hit PMT's. 
The number of hit PMT's within a tight association radius from the projected 
track is counted and the number of photo-electrons measured in the PMT's 
in the association radius is summed. In addition, a quality factor based upon 
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Fig. 17. The electron signal observed in PHENIX. The highest histogram is the 
distribution of the energy and momentum ratio of all charged tracks in 0.8 < pt < 0.9 
GeV/c. The middle histogram is the one with the RICH hits. The lowest histogram 
is the distribution of a random background estimated from the data. The crosses are 
the difference of the middle and the lowest, resulting in the net electron signal. The 
electron signal below E/p < 0.8 are caused by electrons from photon conversions 
far from the beam pipe. 



the position of the associated hits and the position of the projected track is 
calculated. All of these quantities are then used to separate electron tracks 
from the hadron background. 

The gain of the RICH PMT's is determined from the pulse height of the single 
photo-electron peak in the data. The geometry of the mirror system and the 
optics are calibrated from the data using zero magnetic field runs. In these 
runs, charged particle tracks are reconstructed as a straight line connecting the 
event vertex (determined by the BBC), the PCI hit position, and the PC3 (or 
EMCAL) hit positions. The orientation of the mirror panels in the geometry 
class are then adjusted so that the projection point of the track reflected by 
the mirrors points to the center of the Cherenkov rings recorded in the PMT 
array. Fig. 16 shows the distribution of the hit PMT's around the projected 
ray after the geometry calibration. 
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The hadron rejection factor provided by the RICH strongly depends upon 
the multipHcity of the event. From data taken with a test beam with a small 
proto-type of the RICH detector, it was shown that the RICH detector can 
provide a pion rejection factor better than 10^ with an electron efficiency close 
to 100%. However, in high multiplicity events, the electron-pion separation is 
reduced since a hadron track can be associated with a RICH hit produced 
by background electron tracks, resulting in a track that is mis-identified as 
an electron. The pion rejection factor in most central Au+Au collisions was 
determined to be on the order of several hundred with an electron efficiency 
of ~ 80%. Fig. 17 illustrates electron identification using the RICH detector. 
In this figure, the ratio of the energy (measured by the calorimeters) and the 
momentum of the charged tracks in the pt range 0.8 < pt < 0.9 GeV/c are 
plotted. For all charged tracks, the E/p ratio has a mip peak aX, E /p ^ 0.3 
on a broad distribution caused by hadronic interactions. When RICH hits are 
required, a clear peak appears at E/p — 1.0, which is the electron signal. 
The background due to random associations of the RICH and the track is 
estimated, and also shown in Fig. 17. 



10 Time-Of-FIight Detector 



Prior to the application of the Time-Of-Flight detector (TOF) for particle 
identification, the fiight time of a particle impacting the TOF is determined 
by taking the average of the measured stop time values at both ends of the hit 
PMT after subtraction of the start timing provided by the BBC. The pulse 
height information from the PMT's is used to reconstruct the hit position by 
calculating the center-of-gravity of the pulse heights from each end of the slat 
after applying timing calibrations, gain calibrations, geometrical alignment 
corrections, and pulse-height- dependent timing, or slewing effect, corrections. 

For TOF event reconstruction, four distinct C+-|- classes are Tised. The first 
class defines an address object, which provides the mapping between the actual 
read-out channel number of the Front End electronics Module (FEM) and a 
sequential slat identification number. The second class defines a geometry 
object, which handles the conversion from the slat identification number to 
the geometrical hit position of the particle within the PHENIX coordinate 
system. The third class defines a calibration object, which deals with all of 
the calibration parameters used for event reconstruction for the TOF system, 
such as the slat-by-slat timing offset, the position offset, and the PMT gain. 
Finally, the address, geometry, and calibration objects are accessed within a 
TOF reconstruction object within which the energy loss in the scintillators 
and the hit positions are reconstructed. 

During Run 2000, the TOF detector timing offset was cahbrated on a slat- 
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Fig. 18. The timc-of-flight resolution in the transverse momentum range 
1.85 < pt < 1.95 GeV/c for positively charged pions. An overall time-of-flight 
resolution of 115 ps is achieved. 

by-slat basis using reconstructed track information. An overall time-of-flight 
resolution of ~110-120 ps is achieved by selecting pions at high momentum 
and comparing their time-of-flight to that expected from av — c particle with 
a linear trajectory originating at the event vertex, as shown in Fig. 18 after 
corrections for the slewing effect and timing drift during the run. The method 
of particle identification using the TOF detector will be discussed in Section 
13. 



11 Momentum Reconstruction and the Track Model 



Due to the complicated, non-uniform shape of the focusing magnetic field 
along the flight path of charged particles traversing the PHENIX central arm 
spectrometers, an analytic solution for the momentum of the particles cannot 
be determined. Therefore, other approaches such as look-up tables, must be 
used [3]. For Run 2000, a four-dimensional field-integral grid was constructed 
for momentum reconstruction using the drift chamber and for track model def- 
inition within the entire radial extent of the central arms. The variables in the 
field- integral grid are the z coordinate of the event vertex; the polar angle, 
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of the particle at the vertex; the total momentum of the particle, p; and the ra- 
dius, r, at which the field-integral f{p, r, 6q, z) is calculated. The field-integral 
grid is generated by explicitly swimming particles through the measured mag- 
netic field map and numerically integrating to obtain /(p, r, 6'o, z) for each grid 
point. 

An iterative procedure is used to reconstruct the momentum of a reconstructed 
track, utilizing the fact that /(p, r, varies linearly with the angle of the 
track at a given radius. This can be expressed as </) = (l)o + q ■ f{p,r,9o, z)/p. 
Each track is assumed to be a primary track originating from the event vertex 
as determined by the BBC. An initial estimate of the track momentum and 
charge is made from the reconstructed bend angle, a, of the track in the 
drift chamber. The measured polar angle, 6, of the track in the non-bend 
plane at the drift chamber reference radius is used as an initial estimate of 
9q. Then, using the radial position of each reconstructed hit associated to the 
track, a four-dimensional polynomial interpolation of the field-integral grid is 
performed to extract a value of f{p, r, 60, z) for the drift chamber hit. Once this 
is done for all hits, a robust fit in and /(p, r, 6*0, z) is performed to extract 
the quantities 0o and g/p for the track. The extracted values are then fed back 
into the above equation. The initial polar angle, is also determined using 
an iterative procedure using the equation 6 = 60 + 5{p,r = 9o, z) — g{p,r — 
R, 6q, z)/p where 5 is the bend angle of the particle trajectory relative to the 
straight-line trajectory of an infinite-momentum particle. Typically, less than 
four iterations are necessary for convergence on these quantities. 

This procedure has an additional advantage in that it can be used to define the 
shape of the track within the central arm magnetic field, which can be then 
be used to determine the track intersections with each detector in order to 
facilitate inter-detector hit association as described in the next section. This 
is done by storing the coordinates of the particle in radial steps as additional 
entries (but not keys) in the field-integral grid. Line segments connecting the 
interpolated grid coordinates for a track are intersected with the geometry 
objects describing the position of each detector in order to estimate the pro- 
jection of the track on each detector. These projections are used as a seed 
for inter-detector hit association. Finally, the length of the interpolated fine 
segments from the event vertex to a given detector can be summed in order 
to provide an estimate of the fiight distance of the particle to that detector. 
This quantity is used to facilitate particle identification using the TOF. 

The momentum reconstruction resolution is determined by generating sim- 
ulated muons (which will not decay as they traverse the spectrometer) with 

momentum ppe„ and event vertices ranging from —40 < z < 4Q cm. A complete 
GEANT and detailed detector response simulation on each track is performed 
to obtain the detector hits as they would appear in the data, after which 
the reconstruction is run in the same manner as for data, producing a recon- 
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Fig. 19. Momentum resolution as a function of momentum from simulated single 
particle events. 

structed momentum, p^ec- The quantity \Prec~Pgen\/Pgen (in percent) is shown 
in Fig. 19. Multiple scattering dominates the resolution at low momenta, 
while detector resolution determines the resolution at higher momenta. In ad- 
dition, the momentum resolution is not multiplicity-dependent, as observed 
by generating the same plot using central Au-|-Au HIJING events. 



12 Inter-Detector Association 

This section will address the methods and algorithms used for inter-detector 
association. The general object-oriented flowchart used for inter-detector, or 
global, reconstruction is shown in Fig. 20. All stages of the event reconstruc- 
tion query the geometry information for the spectrometer, which is stored as 
Detector Geometry Objects. Geometry interactions (intersections, transforma- 
tions, etc.) are handled either within a specific geometry object type, or within 
a G-|— I- singleton class called the Geometron. The output of detector-specific 
reconstruction are stored in objects containing the hit and track information, 
which are stored in the Data Summary Tapes (DST) created after each recon- 
struction pass. The anchor of inter-detector hit association is the track model 
object, which is discussed further below. Once a track model object is created, 
its methods are used to perform hit association in addition to momentum re- 
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Fig. 20. Illustration of the object-oriented approach to inter-detector association in 
PHENIX. 



construction with the assumption that the particle is a pion. The resulting 
association and kinematics information is stored in distinct objects (and in 
the DSTs), which are then fed into the particle identification procedures. For 
particles that have a probability to not be a pion, the momentum can be 
subsequently reconstructed using the updated particle assumption, as its mo- 
mentum will differ from that of the pion due to energy loss as it traverses the 
spectrometer. 

The track model object is implemented in a flexible manner through a base 
class named PHTrack, which defines the basic functionality that must be pro- 
vided by a generic PHENIX track model class. The momentum reconstruction 
algorithm previously described is implemented as a class that inherits from 
the PHTrack class. The most important application of the track model is to 
determine its intersection with the various detectors in order to perform inter- 
detector track and hit associations. Thus, the data members of the PHTrack 
base class include a list of projection points, vectors, and errors at each de- 
tector. Also included is a list of three-dimensional coordinates approximating 
the shape of the track for use in event displays. It is up to the individual track 
models that inherit from PHTrack to provide the methods for the construction 
and projections of a track prior to associating hits to the track. The basis for 
the construction of a track model object is typically, but not restricted to, a 
drift chamber track. During Run 2000, two distinct track models were utilized. 
The first was an analytical linear track model used for data taken with the 
magnetic field turned off, and the second was the one used for momentum 
reconstruction. 
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The inter-detcctor hit association algorithm operates on each PHTrack object 
within an event individually. Additional inputs for hit association include lists 
of the detector hits (or tracks) which are to be associated to the PHTrack 
object. These "hit lists" are constructed and sorted by spectrometer arm and 
by increasing azimuthal angle using a class developed for this functionality 
in order to simplify and speed up processing. Once constructed, the manip- 
ulation of these "hit list" objects are independent of the detector to which 
the hits or tracks belong. Due to the small number yet wide resolution dis- 
parity of detectors that must be associated, the hit association is performed 
using a three-dimensional road algorithm with each detector associated inde- 
pendently. For a given PHTrack object, its projection point at the detector 
being associated is used to define the center of a window spanning the and 
z coordinate directions. Projection error estimates from the PHTrack object, 
which can be momentum-dependent, are then used to define the widths of 
the roads in and z. These widths can be scaled at run-time via user input 
parameters, typically after examining the residuals of the projections. Once 
the window is defined, it's ranges are fed to the "hit list" object, which re- 
turns either a list of all hits within the window, or the closest hit to the center 
of the window within its boundaries. There is one special case of TEC track 
association considered in this algorithm. For this detector, the a angle of the 
TEC track, along with the a angle of the drift chamber track, are correlated 
with the inverse pt of the track. This information can be used to apply an a- 
difference angle cut between the two detector tracks to reduce the list of TEC 
tracks being considered for association. For the analysis of the data taken in 
Run 2000, the accuracy of the projection determination was good enough that 
the closest hit was used in the initial reconstruction pass while retaining the 
option for specific analyses to further optimize the association hit set in later 
passes on the DST's using the "hit list" methods. 




Fig. 21. The hit association efficiency between the drift chamber and PCI (left), and 
the drift chamber and PCS (right) as a function oipt in central HIJING 200 GeV/A 
Au+Au simulated events. 
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The efficiency of the inter-dctcctor hit association, with the association be- 
ing defined as the hit closest to the center of the projection window, has 
been calculated based upon 200 GeV/A Au+Au central events generated us- 
ing the HIJING event generator [15]. These events were processed through 
the GEANT-based simulation of the detector response for all detectors in the 
central arm spectrometers. The output of the response simulation was then 
reconstructed in the same manner as the data while storing the pointers that 
relate reconstructed global tracks back to the GEANT tracks that contribute 
to them. Utihties that trace a reconstructed object back to the GEANT in- 
formation are implemented as a singleton class for use by a wide variety of 
evaluations. For the evaluation of the global association, the dominant contrib- 
utor of the drift chamber track was used as a basis. The dominant contributor 
is defined as the input GEANT track which contributes the most hits to the 
reconstructed track. In order for an association to be considered to be made 
correctly, the GEANT track that provided the reconstructed hit or track for 
the detector in question must match the dominant contributor of the recon- 
structed drift chamber track to which it is associated. The efficiency is defined 
as the ratio of the number of tracks that were correctly reconstructed and 
geometrically reconstructable to the number of tracks that were geometrically 
reconstructable by each detector (in order to isolate the efficiency to that of 
the association only). The results of the efficiency calculation are shown in 
Fig. 21 for associations to PCI, which is situated adjacent to the drift cham- 
ber, and for associations to PCS after spanning the tracking-detector-free gap 
occupied by the RICH. Although the efficiencies in the left panel of Fig. 21 
are lower, the performance is such that good particle identification can be 
achieved using the TOF, as demonstrated in the next section. 



13 Pcirticle Identification Using the Time-Of-Flight Detector 

Particle identification for charged hadrons is performed by combining the in- 
formation from the drift chamber, PCI, BBC, and the TOF. The time-of-fiight 
resolution of the TOF is measured to be 115 ps, resulting in PID capability for 
high momentum particles. A Aa tt/K separation at momenta up to 2.4 GeV/c, 
and a K/proton separation up to 4.0 GeV/c can be achieved as shown in Fig. 
22. 

Once a global track is reconstructed with information from the drift chamber, 
PCI, and TOF (using a window for TOF association adjusted so that the 
residuals between the projection point and the reconstructed TOF hit position 
is within a standard deviation of 2.5), the flight-path length of the track from 
the event vertex to the TOF as calculated by the momentum reconstruction 
algorithm is used to correct the time-of-flight value measured by the TOF. 
Fig. 22 shows a scatter plot of corrected time-of-flight as a function of the 
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Fig. 22. Scatter plot of the time-of-flight versus reciprocal of the momentum in 
minimum bias Au + Au collisions. This demonstrates the clear particle identification 
capability using the TOF in the Year 2000 data taking period. The flight path is 
corrected assuming the mass for each particle species. 




Fig. 23. The mass squared distribution for positively charged particles integrated 
over all momenta. 



32 



reciprocal of the momentum in minimum-bias Au+Au collisions, after the 
momentum-dependent residual cut between the track projection point and 
TOF hits. The flight path is also corrected for each particle species in this 
figure. Fig. 23 shows the mass-squared distribution for positively charged 
particles integrated over all momenta. The vertical axis in this figure is in 
arbitrary units. It is demonstrated that clear particle identification using the 
TOF is achieved in Run 2000. 



14 Summary 

This document has described the methods and algorithms used for event recon- 
struction in the PHENIX central arm spectrometers. Using a modular, object- 
oriented implementation, widely disparate algorithms from widely disparate 
detector types have been integrated into a coherent reconstruction chain in 
order to successfully reconstruct hits and tracks within each detector, asso- 
ciate hits from different detectors, reconstruct particle momenta, and provide 
particle identification for PHENIX physics analysis. 
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